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Almost all climate models in Coupled Model Inter-comparison Project phase five (CMIP5) were found to 
have a cold bias in Sea Surface Temperature (SST) over the northern Arabian Sea, which is linked to the 
biases in the Indian Summer Monsoon (ISM). This cold SST bias was attributed to the anomalous cold winds 
from the north-western part of south Asian landmass during boreal winter. However, the origin of the 
anomalously strong cold winds over the Arabian Sea and its association with the large-scale circulation is 
obscure. Here we show that an equatorward bias in subtropical Jetstream during boreal spring season 
anomalously cools down the northern Arabian Sea and adjoining land regions in CMIP5 models. The 
models with stronger equatorward bias in subtropical jet are also the ones with stronger cold SST bias over 
the Arabian Sea. The equatorward shift coupled with enhanced strength of the subtropical jet produce a 
stronger upper tropospheric convergence, leading to a subsidence and divergence at lower levels over the 
Arabian deserts. The low entropy air flowing from the Arabian land mass cools the northern Arabian 
Sea. The weaker meridional temperature gradients in the colder models substantially weaken ISM 
precipitation. 

The Sea Surface Temperature (SST) biases are arguably the most prominent error in Coupled General 
Circulation Model (CGCM) simulations 1 , which can result in amplification of model error due to the 
feedback between different components of climate system 2 . Experiments with Atmospheric General 
Circulation Models coupled to prescribed SSTs shows that the external radiative forcing do not directly warm 
up the continents; rather it warms up the oceans which in turn results in a continental warming 1 . This suggests 
that SST biases can have far reaching impact on the simulation of continental climate 4 ' 5 . Thus understanding the 
origins of SST bias is important to improve the model simulations. Tropics-wide bias in SSTs in the fifth phase of 
Coupled Model Inter-comparison Project 6 (CMIP5) models has been traced to biases in the simulations of 
clouds and thermocline depth by the coupled models 7 . Both the local and large-scale oceanic and atmospheric 
processes are dominant elements for the SST biases 8-10 in various ocean basins. Understanding the factors 
responsible for SST biases specific to different ocean basins are vital as those aid climate modelers to rectify these 
biases. In this context, finding the origin of SST biases over the Indian Ocean assumes significance. Further, such 
an analysis also helps in gaining in-depth understanding on the interdependence of various fields in a climate 
model. 

The ocean-atmosphere coupling is very important for the existence of Indian Summer Monsoon (ISM) and to 
maintain its inter-annual and intra-seasonal variability 1112 . Further, precipitation in the monsoon region is 
sensitive to the tropical SSTs 1314 . An accurate simulation of the ISM precipitation is still a challenge to the current 
generation CGCMs 1516 , with the Arabian Sea SST cold bias during the pre-monsoon season highlighted as a major 
reason for the land precipitation error 1719 . The anomalous cold winds from the north-west India and adjoining 
land regions in the CGCMs are suggested as the cause of anomalous cooling of Arabian Sea SSTs 18 . The cold SSTs 
over the Arabian Sea can reduce the meridional temperature gradients that drive the monsoon circulation. 
However, the model mechanisms that cause anomalous cooling of air over south Asian land mass in the pre- 
monsoon season are not yet clear. One possible explanation is that, inadequately resolved orography west of Tibet 
may cause dry air intrusion which is linked to a thermodynamic bias in ISM simulated by many CGCMs 20 . An 
equatorward latitudinal bias in the eddy driven jetstreams in the Southern Hemisphere is identified as a common 
problem in many CGCMs 21 , which is attributed to the biases in midlatitude cloud forcing 22 . Further, asymmetric 
warming between the hemispheres can anomalously shift both Inter-Tropical Convergence Zone (ITCZ) and 
jetstreams towards warmer hemisphere 23 . The double ITCZ, a longstanding problem associated with most of the 
CGCMs, is found to be a response to the increased Southern Hemispheric shortwave forcing 24 . 

Although the equatorward bias in the eddy driven jetstreams and ITCZ are known issues in CGCMS, the 
presence of such biases in the Sub-Tropical Jetstreams (STJ) are not clear. During boreal winter and spring 
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Figure 1 | Arabian Sea SST bias. Northern Arabian Sea (15°N-25°N; 60°E-70°E) spring season (March - May) SST bias (K) calculated for 44 
CMIP5 models. The bias in historical all forcing simulations are computed with reference to fifty year ( 1 95 1-2000) climatology of HadlSST 1.1. This figure 
is plotted using NCAR Command Language (NCL). 



seasons STJ core is located around 30°N which recedes to north of 
—40 °N by the beginning of the summer, in tune with the ISM onset 25 . 
It is possible that the inadequately resolved orography west of the 
Tibetan Plateau in some models 20 may also cause a southward shift of 
STJ over the Southeast Asia. The cold dry air associated with STJ can 
cool the South Asian land mass and the adjacent oceans during pre- 
monsoon season. The anomalous cooling over the ocean can persist 
even after the STJ moves north prior to monsoon onset, owing to the 
higher heat capacity of the ocean. The persisting cold anomaly over 
the northern Indian Ocean can hamper the development of mon- 
soon in its early phase. However, the role of such large-scale circula- 
tion biases in producing anomalously cold SSTs in the Arabian Sea is 
not examined. 

Here we explore the pre-monsoon SST bias over the Arabian Sea in 
CMIP5 historical simulations. We specifically investigate the pos- 
sible presence of an equatorward bias in STJ, in the wake of known 
biases in ITCZ and eddy driven jets 23,24 . 



Results 

All 44 CMIP5 coupled models (see X axis in Fig. 1) analyzed here 
exhibit cold SST bias over the northern Arabian Sea during pre- 
monsoon season (March - May; Fig. 1). The magnitudes of cold 
SST bias in the models range between 0.5 and 3.8 K, consistent with 
the previous analysis 18 . The pre-monsoon climatology of the surface 
temperatures from ERA Interim (ERAI) reanalysis 26 shows Arabian 
Sea SSTs in 299 - 302 K range. A composite of five less cold (WARM; 
Fig. 2b) models show slightly cooler SSTs over the Arabian Sea, 
ranging between 298 and 302K. The composite of five coldest 
(COLD; Fig. 2c) models simulate much colder SSTs. The difference 
between COLD and WARM models clearly shows a strong negative 
gradient in SST over the northern Arabian Sea (Fig. 2d). The contrast 
in surface temperature is even stronger over the land mass adjoining 
the Arabian Sea. The ERAI land surface temperature over the north- 
west India and Arabian deserts exceeds 308K in the pre monsoon 
season (Fig. 2a). Both the WARM and COLD models simulate colder 
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Figure 2 | Pre-monsoon surface temperature and monsoon precipitation. Climatological (MAM) surface temperature (K) from (a) ERAI, (b) WARM 
composite, (c) COLD composites, and (d) difference between COLD and WARM model composites; climatological (JJAS) precipitation (mm day -1 ) 
from (e) GPCP, (f) WARM and (g) COLD model composites and (h) difference between WARM and COLD models. WARM and COLD models are 
indicated in Fig. 1 . GPCP precipitation climatology is calculated for 1 980-2009 period. All other calculations are based onl951-2000 period. This figure is 
plotted using NCL. 
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land surface compared to the observations, with the latter showing 
stronger cooling. The development of ISM circulation is critically 
dependent on the land-sea temperature contrast 11 . The observed 
June - September (JJAS) precipitation (Fig. 2e) shows strong oro- 
graphic precipitation band over the Western Ghats and Arakan 
mountains. The precipitation pattern in WARM models (Fig. 2f) is 
comparable with the observations, although the models simulate less 
precipitation over the Western Ghats. Also, the precipitation band 
extends to northwest India in the WARM models as in observations. 
In the COLD models (Fig. 2g), the JJAS precipitation is substantially 
weaker all over the Southeast Asia compared to WARM composite 
and the precipitation band has not reached northwestern parts of 
India. The difference between COLD and WARM composites 
(Fig. 2h) clearly reveals the weak precipitation in the former, in line 
with previous analysis 15 . This suggests that negative surface temper- 
ature gradients in the COLD models might have slowed down the 
advancement of monsoon circulation. It is worth noting that none of 
the CMIP5 models has been successful in simulating all features of 
ISM satisfactorily 1516 , in line with the cold SST bias. Identifying the 
processes responsible for pre-monsoon cold SST bias may have far 
reaching implications towards an improved model performance of 
ISM simulation. So far, the investigations on cold SST bias are cen- 
tered around local processes, like advection of the cold winds from 
northwest India and Pakistan during boreal winter and spring sea- 
sons 1718 . However, these studies stop short of explaining the origin of 
anomalous cold air over Asian land mass during boreal winter and 
spring. 

A main circulation feature over the subtropics is the STJ, which is a 
response to the equator-pole atmospheric temperature gradient 27 . 
Typically STJ recedes to north of about 40°N during boreal summer, 
in tune with the onset of ISM 25 . The difference in 200 hPa zonal 
winds between COLD and WARM models (Fig. 3a) demonstrate 
that the COLD models have the core of the jet located much closer 
to the equator than WARM models. Over the northern Arabian Sea, 
wind speed at 200 hPa exceeds 8 m s~ ! in COLD models when 
compared to the WARM composite. A comparison with ERAI rea- 
nalysis (Fig. 3b) shows that, the STJ core in COLD models is located 
equatorward (see blue dot in Fig. 3b) of both ERAI and WARM 
models. Further, STJ core is shifted vertically downward toward 
higher pressure level when compared to the ERAI STJ location. In 
contrast, the STJ core in WARM models is located upward of its 
climatological position in ERAI. Further, the strength of STJ core 
is overestimated (slightly underestimated) in COLD (WARM) mod- 
els, with maximum wind speeds of 41 m s -1 (34 m s -1 ) as compared 
to ERAI wind speed of —36 m s~\ The core of STJ is located at 30°N 
at an altitude of 175 hPa in ERAI. The mean position of STJ core in 
WARM (COLD) models is located at 29.5°N and 160 hPa (27.7°N 
and 190 hPa). A latitude-height view of the COLD-WARM zonal 
mean zonal winds (Fig. 3c) clearly shows that the equatorward bias in 
the zonal winds in COLD models is not confined to upper levels, but 
exists throughout the depth of the troposphere. 

The prevalence of strong cold air advection over the Arabian Sea 
can result in anomalous cooling of SSTs in COLD models. In order to 
verify the relation between the anomalous STJ and SST bias over the 
Arabian Sea, a correlation analysis between the climatological (1951- 
2000) zonal winds at 500 hPa and SST bias is carried out for the 
MAM season (Fig. 4a). Both the winds and SSTs are area averaged 
over 15°N-25°N and 60°E-70°E, where the strongest cold anomaly 
is seen over the Arabian Sea. A strong correlation of —0.74 is 
obtained between the two parameters, further strengthening our 
argument that the STJ and pre-monsoon Arabian Sea SST bias are 
related. It is rather straightforward to elucidate the cause-effect rela- 
tionship between the biases in STJ and Arabian Sea SST. Since STJ is a 
planetary scale feature, and the equatorward bias in the COLD mod- 
els span around the globe, it is safe to argue that the SST bias over a 
small region like Arabian Sea may not cause the bias in STJ. On the 



other hand, the analyses presented so far indicate that the bias in 
Arabian Sea SST is most likely linked to the STJ biases. 

The STJ is located on the poleward flanks of the Hadley Cell (HC) 
and the distance between the two STJs on both sides of the equator 
may be considered as a measure of the width of the tropics 28 . The 
composite climatology of the zonal mean meridional mass stream 
function (%) in WARM models shows the descending branch of the 
HC between 5 and 25°N, with its core positioned around mid- 
troposphere (Fig. 4b). On the other hand the HC in the COLD 
composite has a broader and shallow descending structure that spans 
between Equator and 25°N, (Fig. 4c), consistent with the slightly 
downward location of the maxima of STJ (see Fig. 3b). In order to 
bring out the covariability of HC and STJ locations in model simula- 
tions, a linear regression of the zonal mean zonal winds on the pole- 
ward edge latitude of HC (C> H c) across all models is performed 
(Fig. 4d; see methods for details). The poleward latitude of HC is 
taken as the latitude at which % at 500 hPa becomes zero in the 
northern hemisphere. The COLD (WARM) models are found to 
have a mean value of 3> H c = 25.5°N (27.1°N). The latitude-height 
structure of the regression slope indicate that the models in which 
HC subsidence occur at more southward latitude also have the west- 
erly wind maxima located more equatorward. The negative (positive) 
values of regression slopes show linear relationship between 
increased (decreased) westerly winds and more equatorward (pole- 
ward) latitude of HC subsidence. Thus the model-to-model variance 
in the location of STJ maxima is explained by that of Ohc- A ther- 
modynamic bias over the ISM region in climate models was sug- 
gested to be linked to inadequately resolved orography over 
north-western part of India 20 . However, the effect of model resolu- 
tion on the bias in STJ in CMIP5 simulations is not clear. It may be 
noted that MRTCGCM3, which is one of the highest resolution 
(1.125°X 1.125°) models in CMIP5, has the strongest pre-monsoon 
SST bias (see Fig. 1). 

How does the equatorward shift in the STJ produce cold SST bias 
over the Arabian Sea? The mechanism that couples shift in STJ 
location and cooling of the Arabian Sea is explained in Fig. 5. STJ 
causes an upper level positive vorticity, as shown in Fig. 5a, in which 
the MAM climatological composite wind vectors and relative vor- 
ticity at 250 hPa for the WARM models is depicted. The stronger 
winds and its possible interaction with the orography due to a south- 
ward shift of STJ cause stronger relative vorticity in COLD models 
(Fig. 5b). The difference between COLD and WARM models clearly 
shows a stronger upper level convergence over a large region north of 
the Arabian Sea in COLD models (Fig. 5c). The upper level conver- 
gence can induce a subsidence which in turn produces a divergence at 
lower levels. This mechanism is very clear with weak divergence at 
850 hPa in WARM models (Fig. 5d) and strong divergence in COLD 
models (Fig. 5e). The larger values of geopotential height in COLD 
models indicate stronger subsidence at lower levels, as compared to 
WARM models. The low level subsidence location is found to be 
southwest of the upper level convergence area, which may be due 
to the orographic forcing just beneath the upper tropospheric con- 
vergence region. The difference between COLD and WARM com- 
posites shows stronger subsidence and divergence at 850 hPa over 
the Arabian Peninsula and a stronger southward wind flow over the 
Arabian Sea in COLD models (Fig. 5f). This mechanism explains the 
origin of anomalous wind flow over the Arabian Sea in COLD mod- 
els during pre-monsoon period. The wind vectors and equivalent 
potential temperature (8 e ) at 925 hPa show low entropy advection 
towards Arabian Sea in WARM composite (Fig. 6a). The advection of 
low entropy air in COLD models is substantially stronger and 
intrudes further south over the Arabian Sea (Fig. 6b). This low level 
advection of dry and cold air over the Arabian Sea can cool the sea 
surface. Thus the southward shift in the STJ can explain the cooling 
of northern Arabian Sea SST in CMIP5 models. 
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Figure 3 | Horizontal and vertical structure of zonal winds, (a) Difference between COLD and WARM model composites for zonal winds (m s -1 ) at 
200 hPa. (b) Latitude - height cross-section of zonal mean (50°E-100°E) zonal winds from ERA- Interim reanalysis; the red (blue) dots show the mean 
location of jet core in WARM (COLD) models and the vertical and horizontal lines indicate the standard deviation of the core locations. The mean 
latitude ofSTJ core in WARM (COLD) models is 29.5°N ± 2.4° (27.7°N ± 1.8°). The mean altitude ofSTJ core in WARM (COLD) models is 160 ± 22.4 
(190 ± 22.4) hPa. The mean ± standard deviation of jet core speed in WARM (COLD) models is 34 ± 3.5 (41 ± 3.9) m s _l . The black dot represents the 
core of jet in ERAI (30°N and 175 hPa). The mean jet core speed of ERAI is 36 m s _1 ; (c) Latitude - height cross-section of COLD - WARM model 
composites of zonal mean zonal winds. All calculations are done using fifty years (1951-2000) climatological mean spring season, except in the case of 
ERAI, for which thirty year (1980-2009) climatology is used. This figure is plotted using NCL. 



Discussion 

We have presented a mechanism that explains the Arabian Sea cold 
SST bias during pre-monsoon season in many CGCMs. The current 
understanding on the issue tries to attribute the Arabian Sea cold SST 
bias to the anomalous cold air advection from northwest India and 
adjoining land regions. However, these investigations fall short of 
identifying the source of anomalous cooling of continental air over 
the Asian land mass. It is important to understand the exact cause of 
the SST bias in order to improve the performance of the CGCMS. 



Our analysis finds that the origin of anomalously cold air over the 
South Asian land mass in CGCMs originates from a southward shift 
coupled with increased strength of the STJ. The biases in the location 
and strength of STJ are explained by the location of northern hemi- 
spheric Hadley Cell subsidence. The biases in the strength and loca- 
tion of Hadley Cell subsidence may be linked to the biases in the 
radiative forcing 24 . Our results demonstrate the linkage between the 
large-scale circulation and a local phenomenon at a relatively smaller 
ocean basin and the resultant impact on regional climate. The biases 
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Figure 4 | Wind speed versus SST bias, (a) Scatter plot between area averaged (15°N-25°N and 60°E-70°E) 500 hPa zonal wind (m s" 1 ) and 
northern Arabian Sea SST bias (K) during March - May Season; Hadley Cell and zonal winds, (b) WARM and (c) COLD model composites of zonal mean 
stream function (10 10 kgs" 1 ); (d) linear regression of Hadley Cell latitude on zonal mean (50°E-100°E) zonal winds during March - May season. Stippling 
in (d) shows statistically significant (p<0.05) regression slopes as revealed by a two-tailed t-test. This figure is plotted using NCL. 



in SST and STJ are found to be independent of the temporal scale 
used for the calculations (see methods for details). Since the summer 
monsoon precipitation influences the life of the entire South Asian 
continent with a massive population, prediction of monsoon rainfall 
assumes significance. The results in this study help to identify the 
cause of the most prominent systematic bias seen in climate models 
in simulating the seasonal mean monsoon precipitation. 

Methods 

We used monthly mean outputs of CMIP5 historical all forcing simulations. The 
period during 1951-2000 is used for the calculations presented here. The fifty year 
climatology is sufficient to address the natural variability of various temporal fre- 
quencies. Our focus is pre-monsoon season, as the strong ocean-atmosphere coupling 
during monsoon season makes it difficult to attribute the SST bias to a source other 
than the coupled processes. The Hadley Center Sea Surface Temperature 29 
(HadISSTl.1) is used to calculate the SST bias in CMIP5 models. The observation 
density in HadISSTl.1 in the second half of the 20 th century is better compared to the 
first half 29 , which is the reason we choose 1951-2000 as the base period for our 
calculations. However, in order to ensure robustness, the calculations are repeated for 
another two temporal scales (not shown here) - entire 20 Lh century and 1981-2000 
period. All the three temporal scales chosen for the analysis yielded identical results, 
suggesting that the biases are independent of temporal scale. The surface temperature 
form model output is used as a proxy for model SST. The SST bias is computed over 
the Arabian Sea region of 15°N-25°N and 60°E-70°E, The zonal mean of the zonal 
winds are computed over the longitudes 50°E-100°E. The zonal winds from 
ERA-Interim (ERAI) reanalysis 26 are used for qualitative comparison of the model 



simulated zonal winds. Also the surface temperature from ERAI is used to compare 
with the model simulated surface temperature. The Global Precipitation 
Climatology 30 version 2.2 (GPCP) monthly precipitation data is used to compare the 
model simulated precipitation. All the model data are interpolated to a common 
2.5°X2.5 LI latitude-longitude grid using linear interpolation, before doing 
calculations. 

The core of subtropical Jetstream is calculated in a way identical to the technique 
proposed to calculate the core of mid-latitude Jetstream 11 . The latitudinal location of 



the core of the jet, <j)- el 



S u>uJU(<j>)d<l>' 



where 0 is the latitude, U x is a threshold 



value of the zonal mean (50°E- 100 E) zonal winds ( U) around the jet core, over which 
the integral is calculated. Here U x is defined as U x — 0.9U max> where U max is the 
maximum zonal wind at the jet core (grid point maximum). The structure of the jet is 
symmetric closer to the core and becomes asymmetric away from it (see Fig. 3b). The 
value of U x should be taken in such a way that the integral is calculated over a 
symmetric curve around U max . Otherwise this method may produce a bias in the 
calculation of the &j et . The advantage of calculating &j et using this method rather than 
taking the latitude of u max is that we can overcome the issue of coarse grid resolution. 
The regression slope in Fig. 4d is calculated as follows. Regression coefficient at a 

particular grid point, /J = ^— UtUjj ^— Ui^A ; where U denotes zonal 

wind at a grid point, 0 the poleward latitudinal boundary of northern hemispheric 
Hadley Cell, the subscript i stands for the index of CMIP5 models (span from 1 to 43) 
and V for total number of models used in the calculation (in this case n— 43). It may 
be noted that the three dimensional monthly mean wind data are not available for the 
model CESM1 -CAM5- 1-FV2 (historical experiment) in CMIP5 archive while car- 
rying out this study, and hence omitted in the regression analysis. The models are 
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Figure 5 | Upper level convergence and lower level divergence. Wind vector (ms _1 ) and relative vorticity (10~ 5 s~') at 250 hPafor (a) WARM, (b) COLD 
model composites, and (c) difference between COLD and WARM model composites; wind vector and geopotential height (m) at850 hPa for (d) WARM, 
(e) COLD model composites, and (f) difference between COLD and WARM model composites. This figure is plotted using NCL. 

organized in alphabetical order of their names. The poleward latitude of HC is the the latitude of Hadley Cell subsidence. Two tailed t-test is used to identify the values of 
latitude at which x at 500 hPa becomes zero. This calculation is carried out for each regression slopes that are statistically significant at 5% level (P-value<0.05). 
latitude-height grid cell to generate the regression map. This regression map explains Graphics software: All plots are produced with NCAR Command 
the relation between model-to-model variability in the zonal mean zonal winds and Language. 
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